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van der Waals coefficients for positronium interactions with atoms 
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The random-phase approximation with exchange (RPAE) is used with a B-spline basis to compute 
dynamic dipole polarizabilities of noble-gas atoms and several other closed-shell atoms (Be, Mg, 
Ca, Zn, Sr, Cd, and Ba). From these, values of the van der Waals Ce constants for positronium 
interactions with these atoms are determined and compared with existing data. Our best predictions 
of Ce for Ps-noble-gas pairs are expected to be accurate to within 1%, and to within a few per cent 
for the alkaline earths. We also used accurate dynamic dipole polarizabilities from the literature to 
compute the Ce coefficients for the alkali-metal atoms. Implications of increased Ce values for Ps 
scattering from more polarizable atoms are discussed. 
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I. INTRODUCTION 

The interaction of positronium (Ps) with matter and 
antimatter is an important topic [1] with applications in 
many areas of physics. For instance, the proposed AEgIS 
experiment would produce an antihydrogen beam from 
the reaction between Ps and antiprotons [2, 3]. The anti¬ 
hydrogen would then be used to determine whether an¬ 
timatter is affected by gravity in the same manner as 
matter. Ps is widely used in condensed matter physics to 
determine pore sizes in microporous materials and probe 
intermolecular voids in polymers [4]. Further, positron¬ 
ium formation in porous materials is used to study its 
interactions with gases, e.g., xenon [5, 6], or the interac¬ 
tion between the Ps atoms themselves, with prospects of 
Bose-Einstein condensation at room temperature [7, 8]. 

Here we focus on the problem of Ps-atom scattering. 
Compared with electron-atom scattering and positron- 
atom scattering, Ps-atom scattering is more difficult to 
treat theoretically, chiefly because both scattering ob¬ 
jects have an internal structure [9]. 

In this paper we address low-energy Ps scattering from 
closed-shell atoms. The short-range Ps-atom interaction 
is repulsive, because (a) the Pauli principle prevents the 
electron from Ps from entering the volume occupied by 
the atomic electrons, and (b) the positron is repelled by 
the screened potential of the atomic nucleus. However, 
low-energy Ps-atom scattering is also affected by the at¬ 
tractive long-range van der Waals interaction [10-12]. 

The van der Waals potential behaves asymptotically 
as 

U(R) ~ -^§, ( 1 ) 

where R is the distance between the centers of mass of 
the atoms, and Ce is the van der Waals coefficient for the 
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atomic pair [13]. The value of the Cq constant determines 
the scattering phase shifts for the partial waves l > 2 at 
low collision momenta k [14], 


6C 6 7rfc 4 

U ~ (21 - 3)(2Z - 1)(2Z + 1)(2Z + 3)(2 1 + 5)' 


( 2 ) 


The magnitude of Cq also affects the Ps-atom scattering 
length A. This can be seen from the estimate which uses 
the potential (1) with a cut-off at R = Rq [12], 

/mG 6 \ 1/4 r(3/4) J_i/ 4 (x 0 ) 
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where m is the reduced mass, T and J„ are the gamma 
and Bessel functions, respectively, and 
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( 4 ) 


The dimensional prefactor in Eq. (3) determines the char¬ 
acteristic magnitude of the scattering length in atomic 
collisions [15]. We use atomic units throughout, so that 
m = 2 for Ps-atom collisions. 

Mitroy and Bromley [16] calculated the Ce constants 
for Ps-noble-gas interactions using the oscillator strength 
sum rule. For He they calculated the oscillator strengths 
via the configuration interaction (Cl) method, while for 
Ne, Ar, Kr, and Xe recourse was made to a set of pub¬ 
lished pseudo-excitation energies and dipole oscillator 
strengths [17], leading to semi-empirical (SE) values of 
Ce- In a recent paper devoted to Ps scattering from Al¬ 
and Kr [12], Fabrikant and Gribakin used the London for¬ 
mula [18] for the van der Waals coefficients, which gave 
values about 6% greater than those from Ref. [16]. To es¬ 
timate the effect of this difference on Ps-atom scattering 
we can use Eq. (3), which shows that for Ar (i?o = 2.67 
a.u. [12]), changing Ce from 152 a.u. [12] to 144.1 a.u. [16] 
increases the scattering length by about 4%. This exam¬ 
ple shows that Ps-atom scattering is sensitive to the value 
of the van der Waals coefficient. 

In this work, we employ the random-phase approxima¬ 
tion with exchange (RPAE) [19, 20] to calculate the dy¬ 
namic dipole polarizabilities of the noble-gas atoms and 
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several other closed-shell atoms with an ns 2 valence shell. 
These are used to calculate the C e constants for Ps-atom 
pairs ab initio. The RPAE method, often called simply 
the random-phase approximation (RPA), is equivalent 
to the (linearized) time-dependent coupled Hartree-Fock 
method (see, e.g., Ref. [21] and references therein). RPAE 
is known to provide a good description of atomic polariza¬ 
tion for the noble gases [22, 23]. To make our predictions 
more accurate, we use scaling based on the known values 
of the static dipole polarzabilities. As a result, we expect 
our final C e constants to be accurate to better than 1%. 
As far as we are aware, there are currently no published 
values of Ce for Ps interactions with other closed-shell 
atoms. Although the RPAE method is less accurate for 
them, our final Ce values provide a good benchmark for 
future work. 


atoms (such as the noble gases) the single-particle (e.g., 
Hartree-Fock) method does not give accurate results. 
RPAE is a many-body theory method that is known to 
give accurate dipole polarizabilities and photoionization 
cross sections for closed-shell atoms, with the best results 
for the noble gases [19, 22], 

After the RPAE equations have been solved (see Ap¬ 
pendix A for details), the dynamic dipole polarizability 
is calculated as 
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II. THEORY 

The van der Waals coefficient for a pair of atoms, A 
and B, may be expressed as 

^ 2^ |(n^|D|0 A )| 2 |(n^|D|0B>| 2 _ 

° 6 3^ E*+E$-E#-E? (5j 

where D is the electric dipole operator, the matrix ele¬ 
ments are taken between the ground (0) and excited (n or 
n') states of the atoms, with respective energies E 0 and 
E n , and the extra indices (A or B) are used to distinguish 
the two atoms [24]. 

There exists a useful relationship connecting dynamic 
polarizabilities of imaginary frequencies with Ce- The dy¬ 
namic polarizability at frequency u> is given by 

a(u) = I V ( IHD|0)l 2 , |(n|D|0)| 2 \ 

3 ^ l E n — Eq — oj — id E n — Eq + u> + iS / 

( 6 ) 

where S is a positive infinitesimal, which determines the 
sign of the imaginary part of a(u) for real frequencies ui 
above the ionization potential of the atom. The ground- 
state static polarizability is a(0). The van der Waals con¬ 
stant can be obtained by integrating the product of dy¬ 
namic polarizabilities of atoms A and B over imaginary 
frequencies, viz. 

3 f°° 

Ce = — a.A(ioj)aB{iu) dio. (7) 

7T JO 

This result is very convenient because the integration 
path avoids the poles of the dynamic polarizabilities. 

For hydrogenic atoms (such as H and Ps), the polar¬ 
izabilities are given by Eq. (6) with the dipole matrix 
elements replaced by the single-particle radial matrix el¬ 
ements (np|r|ls), where r is either the electron-proton 
(H) or electron-positron (Ps) separation. Thus, for an in¬ 
teracting pair of hydrogenic atoms, an essentially exact 
value of Ce can be found. However, for many-electron 


where (i/||d||/i) is the reduced Hartree-Fock dipole matrix 
element, (^||A(u;)||/r) is the reduced RPAE dipole matrix 
element, e„ is the energy of state v, and id has been 
dropped for ui = 0 and imaginary frequencies. In Eq. (8) 
the sums over the magnetic quantum numbers and spins 
have already been carried out, and the remaining sums 
are over the occupied (< F) or empty (> F) electron 
orbitals v and [i. 

Note that the matrix element (^||A(w)||/j,) is calculated 
off-mass-shell, i.e., for ui ^ e„ — e M , and the energy differ¬ 
ences in the denominator of Eq. (8) contain the unper¬ 
turbed Hartree-Fock energies. An equivalent but more 
complicated approach would be to calculate these ampli¬ 
tudes on the mass-shell, simultaneously with finding the 
RPAE excitation energies In this case the expression 
for the polarizability would contain modulus squared val¬ 
ues of (n\\A(oj v ^)\\n) and the RPAE excitation energies 
in the denominator [cf. Eq. (6)]. 

The expressions for Ce and a(u>) contain sums over the 
complete sets of excited states. In real systems, such as 
atoms, these sets of excited states include both discrete, 
Rydberg states and the continuum of states with energies 
above the ionization potential of the system. By using B- 
splines in a box of finite radius R (see Section III), the 
continuum is discretized in a way that allows accurate 
numerical calculations of both the static polarizability 
(for ui = 0) and a(iw), and Ce- 

III. NUMERICAL RESULTS 

We use B-splines to construct either hydrogenic or 
Hartree-Fock basis states. They provide an accurate rep¬ 
resentation of the ground-state orbitals and an effective 
spanning of the continuum, due to an appropriately cho¬ 
sen radial grid (see, e.g., Ref. [25]). In this work we use a 
set of 60 R-splines of order 9, with a box size of R = 15 
a.u. The absence of true continuum states does not affect 
the accuracy of a(0) or a(iui) because all excitations in 
the sums are virtual (the denominator never vanishes); 
such virtually excited electrons cannot travel away from 
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TABLE I. Static dipole polarizabilities a( 0) for the noble 
gases and other closed-shell atoms (in atomic units). 


Atom 

Present 

RPAE a 

RRPA b 

Ree. c 

F d 

He 

1.322 

1.256 

1.322 

1.384 

1.047 

Ne 

2.377 

2.30 

2.38 

2.67 

1.123 

Ar 

10.758 

10.73 

10.77 

11.07 

1.029 

Kr 

16.476 

16.18 

16.47 

17.075 

1.036 

Xe 

27.099 

27.98 

26.97 

27.815 

1.026 

Be 

45.604 

43.2 

45.6 

37.76 

0.828 

Mg 

81.502 

— 

81.2 

71.3 

0.875 

Ca 

183.965 

166 

182.8 

157.1 

0.854 

Zn 

54.046 

— 

50.8 

38.8 

0.718 

Sr 

242.240 

— 

232.6 

197.2 

0.814 

Cd 

75.958 

— 

63.7 

49.65 

0.654 

Ba 

355.735 

— 

324.0 

273.5 

0.769 


a This implementation [22, 26] suffers from poor convergence. 
b Relativistic random-phase approximation [27]. 
c Recommended experimental values for the noble gases, Zn, and 
Cd [28], and calculated values for alkaline earth metals [29]. 
d F is the ratio of the recommended value to the present value. 


the atom, so describing them by a set of states in a box 
is accurate. 

We calculated static dipole polarizabilities of H, Ps, the 
noble gases, and other closed-shell atoms and compared 
with exact theoretical or best experimental or theoreti¬ 
cal values to verify the validity of the method. Then we 
computed dynamic polarizabilties a(iu>j ) over a discrete 
set of imaginary frequencies 

uij =uj 0 [e aij ~ 1) -1], j = l,...,N ul , (9) 

where wo, N u , and cd max = Wj\r w are parameters, and 


a = 


1 

N u - 1 


In 



( 10 ) 


Values of = 0.01 a.u., N u = 100, and w max = 1000 a.u. 
have been used throughout; these values were chosen to 
provide accurate, converged values of the integral (7). For 
H and Ps the polarizabilities were calculated essentially 
exactly; for the noble gases and ns 2 atoms they were 
calculated using the RPAE method. The values of Cq 
were then found by evaluating Eq. (7) numerically. 

The static polarizability obtained for H was a(0) = 
4.500 a.u., in perfect agreement with the exact value of 
a(0) = 9/2 [24], The U-spline states and dipole ampli¬ 
tudes for hydrogen can be used to calculate a(0) and 
a(ico) for Ps by halving the energies and doubling the 
amplitudes (due to the reduced mass of Ps being a half 
of that of H). This gives a(0) = 36 a.u. for Ps. 

Table I shows the static polarizabilities obtained us¬ 
ing RPAE for the noble-gas and other closed-shell atoms. 
The RPAE static polarizabilities for the noble-gas atoms 
agree to 0.1% or better with the results of the equiv¬ 
alent coupled Hartree-Fock calculation [30], with the 
present values for Kr and Xe being more accurate nu¬ 
merically due to better convergence. The results are also 


TABLE II. van der Waals Ce coefficients for various Ps-A' 
systems (in atomic units). 


System 

Present 

Best prediction 3 

CI/SE b 

Other 0 

Ps-H 

34.785 

— 

34.785 

34.785 

Ps-Ps 

207.969 

— 

— 

207.966 

Ps-He 

12.849 

13.41 

13.37 

14.6 

Ps-Ne 

23.759 

26.48 

26.74 

27.4 

Ps-Ar 

96.212 

98.69 

98.50 

104.4 

Ps-Kr 

142.185 

146.71 

144.1 

155.1 

Ps-Xe 

222.355 

227.38 

221.6 

240.6 

Ps-Be 

241.416 

208.6 

210.7 

294.3 

Ps-Mg 

393.484 

355.8 

357.9 

506.4 

Ps-Ca 

717.138 

641.5 

636.6 

1079 

Ps-Zn 

300.756 

231.4 

— 

303.8 

Ps-Sr 

885.425 

767.5 

763.8 

1144 

Ps-Cd 

397.161 

285.0 

— 

381.6 

Ps-Ba 

1164.429 

975.6 

966.8 

1569 


a Obtained by scaling the matrix elements by £ and energies by 
l/e 2 , where £ = F 1 / 4 (see Table I). 
b Computed in Ref. [16] using the configuration-interaction (Cl) 
for He and semi-empirical dipole polarizabilities [17] for the 
noble gases, and in the present work, using empirically-adjusted 
Cl + many-body theory data [31] for the alkaline-earth atoms. 
c For Ps-H and Ps-Ps: pseudostate calculations [32], For all other 
systems, these values were calculated using the London formula 
C& = 3aAas/A^s/[2(/A + 7 b)], where oa,b and Ia,b are the 
static dipole polarizabilities and ionization potentials of the 
atoms [18]; the polarizabilities were taken from Ref. [28] and 
the ionization potentials were taken from Ref. [33]. 


in agreement with calculations made using the relativis¬ 
tic random-phase approximation (RRPA, which also ac¬ 
counts for exchange). Larger differences are to be ex¬ 
pected for heavier atoms, since relativistic corrections 
scale as (Z/137) 2 , where Z is the nuclear charge. How¬ 
ever, even for the heaviest noble-gas atom (Xe) the dif¬ 
ference is less than 0.5%. This bodes well for using the 
dynamic polarizabilities to calculate the Cq coefficients 
for Ps-noble-gas pairs. As expected, for the ns 2 atoms, 
the agreement with RRPA is poorer for heavier atoms 
(since the s electrons are affected more strongly by the 
relativistic corrections). The agreement with the experi¬ 
mental values for these atoms is also poorer in general. 
This is related to the smaller ionization potentials and 
larger effect of the non-RPA correlation effects (e.g., two- 
hole-two-particlc excitations) in these systems. 

The results for Cq are displayed in Table II. Our calcu¬ 
lations which use RPAE polarizabilities (“Present” in Ta¬ 
ble II) are in close agreement with the SE calculations of 
Mitroy and Bromley [16] for the heavier noble gases; the 
relative differences for Ps-Ar, Ps-Kr, and Ps-Xe are 2%, 
1%, and 0.3% respectively. For the lighter noble gases the 
differences are more significant: 4% and 11% for Ps-He 
and Ps-Ne, respectively. This discrepancy can be traced 
back to the fact that the RPAE polarizabilities for He 
and Ne are lower than the recommended values, by 5% 
and 12%, respectively. Correcting the RPAE C$ value 
for He by the corresponding factor F (see Table I) gives 
Cq = 13.45, which agrees with the accurate value for He 
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from Ref. [16] to within 0.6%. However, this crude scal¬ 
ing violates the Thomas-Reiche-Kuhn sum rule, which in 
the limit of large oj gives 

N 

a(iuj) ~ -g, (11) 

or 

where N is the number of atomic electrons. A better ap¬ 
proach is to scale the dipole matrix elements (/j||d||z') and 
(y\\A(u))\\v) in Eq. (8) by some factor £ and the orbital 
energies in the denominators by l/£ 2 . In this case the 
sum rule is preserved, while the calculated static polar¬ 
izability is scaled by £ 4 . The appropriate choice of £ is 
given by £ = F 1 / 4 . Carrying out this scaling for He gives 
Cq = 13.41, which agrees with the value from Ref. [16] to 
within 0.3%. This is a clear improvement over the simple 
scaling by a factor of F. 

As a further test of the effectiveness of this scaling, 
the Cq constants were calculated for pairs of the noble- 
gas and alkaline-earth-metal atoms. These were com¬ 
pared with data from Ref. [31], wherein several relativis¬ 
tic many-body theory [29, 34] and semi-empirical [35] 
methods were employed. For noble-gas pairs, the mini¬ 
mum relative difference was 0.7% (for Ar-Ar), while the 
maximum (for Xe-Xe) was 5% (due to the use of the dy¬ 
namical polarizability “normalized” to a(0) = 27.16 a.u. 
[35] in Ref. [31], which is lower than the recommended 
experimental value a(0) = 27.815 a.u. that we use). For 
alkaline-earth-metal pairs, the comparison was actually 
better; the relative differences ranged from 0.2% (for Sr- 
Sr) to 2.1% (for Be-Be). This is because we use the same 
recommended static polarizabilities for the alkaline earth 
atoms as in Ref. [31]. 

Scaling the ab initio RPAE polarizabilities in this way 
produces our best prediction of the Cq van der Waals 
coefficients (third column in Table II). We expect that for 
the noble-gas atoms these values are accurate to within 
1%. The relative differences from the CI/SE values of 
Mitroy and Bromley for Ne, Ar, Kr, and Xe are 1%, 0.2%, 
1.8%, and 2.6%, respectively. We believe that for Kr and 
Xe our best prediction values are superior to those of Ref. 
[16], where the polarizabilities from Ref. [35] were used. 

Since the RPAE polarizabilities of the ns 2 atoms are 
greater than the recommended values of a(0) by 15-50%, 
the calculated ab initio Cq values significantly overesti¬ 
mate the true van der Waals coefficients. Given the larger 
discrepancy for a(0), the use of scaling is a cruder pro¬ 
cedure for improving the Cq constants. In this case we 
believe that our best predictions are accurate within few 
per cent for the alkaline earth atoms, and within 5-10% 
for Zn and Cd. 

As a further test of the accuracy of our predictions 
for the alkaline-earth atoms, we used tabulated dynamic 
dipole polarizabilities from Ref. [31], which were com¬ 
puted using a combination of relativistic methods, in¬ 
cluding RRPA, Cl and many-body perturbation theory, 
and further adjusted using accurate theoretical and ex¬ 
perimental data. These values are presented in the fourth 
column of Table II (lower half). In all cases the differ- 


TABLE III. Van der Waals Cq coefficients for Ps-alkali-metal 
systems (in atomic units). 


Atom 

Li 

Na 

K 

Rb 

Cs 

Fr 

C 6 a 

C 6 h 

462.8 

463.9 

505.7 

507.2 

765.7 

768.3 

846.8 

844.8 

1014.6 

937.5 


a Computed in the present work using empirically-adjusted 
many-body theory data [31]. 
b Semi-empirical calculations from Ref. [16]. 


ence between these values and our best predictions does 
not exceed 1%. Ref. [31] also presents accurate dynamic 
dipole polarizabilities for the alkali-metal atoms, and we 
used these to compute the Cq coefficients for the Ps- 
alkali-metal pairs, see Table III. Comparison with semi- 
empirical calculations presented in Ref. [16] shows excel¬ 
lent agreement. 

Looking at the last column in Table II, we see that the 
London formula does a reasonable job for the more rigid 
noble-gas atoms, but tends to overestimate the Cq coef¬ 
ficients significantly for Ps interacting with ns 2 atoms. 


IV. CONCLUSIONS 

Dynamic dipole polarizabilities over a range of imagi¬ 
nary frequencies were computed numerically exactly for 
H and Ps, and by using RPAE for the noble-gas atoms 
and several other closed-shell atoms. The static polar¬ 
izabilities for the noble gases deviated from relativistic 
(RRPA) calculations by no more than 0.5%, and for Ar, 
Kr, and Xe were within few per cent of experimental val¬ 
ues. There was greater error for the other closed-shell 
atoms, but this was expected; the RPAE method is most 
suitable for the noble gases. 

Ab initio calculations of the van der Waals Cq coef¬ 
ficients for Ps interactions with these atoms were per¬ 
formed. For the heavier noble gases, close agreement is 
observed with previous calculations [16]. For the lighter 
noble gases (He and Ne) the discrepancies are more sig¬ 
nificant, which motivates a sum-rule preserving scaling 
of the dynamical polarizabilities to calculate Cq values. 
Though these data are no longer ab initio , they are ex¬ 
pected to the be the most accurate values currently avail¬ 
able. For the most part, our values of Cq for Ps inter¬ 
actions with the ns 2 atoms differ significantly from the 
London values. Here our best predicitons are less accu¬ 
rate (few to 10%), but as there are no other theoretical 
calculations of these data, these values will provide a use¬ 
ful benchmark for future calculations. 

It is interesting to discuss the implications of the Cq 
values for low-energy Ps-atom scattering. It is clear from 
Table II that for more polarizable (and more weakly 
bound) atoms, the Cq values are greater, and that Ps 
will experience a stronger van der Waals attraction to 
these atoms. However, it can be seen from the parame¬ 
ter xq, Eq. (4), which determines the scattering length 
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in Eq. (3), that the increase in Cq for such atoms is off¬ 
set by the increase in the parameter R$, which is pro¬ 
portional to the atomic radius. In fact, the latter effect 
makes xo smaller for the more weakly bound atoms. As 
a result, the more polarizable atoms are not more attrac¬ 
tive for Ps, and will likely have larger positive Ps scatter¬ 
ing lengths, due to their larger geometric sizes. This also 
shows that Ps binding to closed-shell neutral atoms does 
not occur. Note that this is in contrast to many open-shell 
atoms, e.g., Na, Cu, or the halogens, which do bind Ps by 
accommodating the extra electron in their valence shell 
[36, 37]. However, even in the case of Ps bound states 
with alkali-metal atoms, the binding energy decreases in 
the sequence PsLi, PsNa, PsK [38], in spite of the greater 
values of the C b constant. This is primarily the effect of 
the increasing atomic radius. 

It is hoped that the results presented here will be use¬ 
ful for studies of Ps-noble-gas-atom scattering and Ps 
interactions with ns 2 and alkali-metal atoms. With lit¬ 


tle extra work, the method can be extended to calcu¬ 
late quadrupole and higher polarizabilities and determine 
higher-order van der Waals coefficients Cg and C\q. 


ACKNOWLEDGMENTS 

The work of A.R.S. has been supported by the Depart¬ 
ment of Employment and Learning (Northern Ireland). 


Appendix A: RPAE equations 

In the equations for the RPAE dipole matrix element 
{v\A{<jj)\n), we must distinguish between the hole states, 
i.e., states below the Fermi level F, and particle states, 
i.e., states above the Fermi level F. For v > F and p < F 
we have 


{v\A(uj)\h) = {v\d\p) 


E 

i y>F,^j'<F 





(V W) - (VI vWn) 

co — e v > + + i( 1 — 2 n v >)6 




(Al) 


and a formally identical equation for v < F and p > F. 

In Eq. (Al), (y\d\p) is the Hartree-Fock dipole matrix element; (vp'\V\v' p) is the Coulomb matrix element, which 
is defined by 

{vp'\V\v'p) = JJ y£(r)p*,(r , ) | r ^ r ,| yv(r')(r) d 3 r ri V, (A2) 

where the ip v are single-particle wave functions; e„/ is the energy of state v'\ and 


0 for v' > F, 
1 for v' < F. 


(A3) 


By separating the angular and radial parts in the electronic states ip, and then integrating over the angular variables 
and summing over the magnetic quantum numbers and spins, one obtains the RPAE equations for the reduced 
amplitudes (^||A(w)||/i) in the form 


where 


(HIAMIIm) 


Hd\\n) + l 


E 

v'>F,n' <F 



Il'>F,v' <F J 


(np'WU.y^WWAjioW) ' 

co — e u ' + + i(l — 2 n v ')5 ’ 


(y\\d\\p) 



P v {r)rP^(r) dr, 


(A4) 


(AS) 


(VIIVIV) = 2(WillV) -sf^-i)*- 1 {/" } M (VINIV), (A6) 

i=o ^ M v > 

and 

l 0 0 J* J o Rp »( r ) P Ar')^PAr')P l >(r)drdr\ (A7) 

are the reduced dipole and Coulomb matrix elements, P„{r) are radial wave functions, [/„] = 2 l v + 1, r> = max(r, r'), 
r< = min(r, r'), and R is the box radius in our R-spline-basis implementation. 



Introducing vectors x and y for (z/||A(w)||p,) for v > F, 
p < F and v < F, p > F, respectively, we can write 


equations (A4) in block matrix form as 




fu la 

u lb \ 

(XI 

0 ] 

\u lb 

u la 

Vo 

X2j 



(A8) 
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where xi and X 2 are the diagonal matrices of energy de¬ 
nominators, d is the vector of Hartree-Fock dipole matrix 
elements, and the matrices U\ a and Uu represent the two 


terms in (A6). This linear equation can be solved nu¬ 
merically for the RPAE dipole amplitudes, and then the 
dynamic dipole polarizability is calculated from Eq. (8). 
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